The normal state scattering rate in high-T c cuprates 
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I present a new phenomenological model for the normal state transport properties of high-T c 
cuprates. In particular, I identify a form of scattering rate that may account, qualitatively 
and quantitatively, for the normal state (magneto)-transport properties of Tl2Ba2CuOe+« and 
Bi2Sr2CaCu208+i from optimal doping through to the overdoped side of the phase diagram. The 
form of the scattering rate is also consistent with features seen in photoemission spectroscopy in 
Bi2Sr2CaCu208+i and offers a new intuitive way to understand the evolution of the temperature 
dependence of the inverse Hall angle with disorder and with carrier concentration. 



PACS numbers: 72.10.-d, 74.25.Fy, 74.72.Fq 
I. INTRODUCTION 

The search for a consistent description for the 
anomalous normal state transport properties of high-T c 
cuprates (HTC) has been a long and fascinating jour- 
ney. Many phenomenological models seeking to identify 
the appropriate form for the normal state scattering rate 
have been offered but none have captured the complete 
picture. From the very beginning, it was clear that the T- 
linear behavior of the in-plane resistivity p a b, extending 
over an anomalously broad range of temperature for op- 
timally doped (OP) cuprates, would be difficult to recon- 
cile within a conventional Fermi-liquid (FL) framework.EJ 
Later observations of the T 2 dependence and, impurity 
dependence of the inverse Hall angle cot©#p together 
with the modified Kohler's rule (in-plane orbital mag- 
netoresistance (MR) Ap a b/p a b ~ tan 2 0#),0 seemed to 
constitute an insurmountable challenge to FL theory. 

Several notable attempts to explain the anomalous 
transport behavior of HTC within a modified FL scenario 
followed, each one incorporating a particular form for 
the (strong) anisotropy of the transport scattering rate 
within the plane arising out of some form of couoiing to a 
singular bosonic mode; be it spin fluctuations,™! charge 
fluctuationsQ or ei-wave superconducting fluctuations.Q In 
each case, however, inconsistencies with certain aspects 
of the experimental data were exposed. 

Given these failings, other more exotic models, based 
on a non-FL ground state, gained prominence within 
the community; most notably the two-lifetime picture of 
AndersonO and the Margiiial-Fermi-liquid (MFL) model 
of Varma and co-workers The two-lifetime model has 
been extremely insightful for understanding the experi- 
mental situation in OP cuprates, but has yet to provide 
an explanation for the results of Angle-Resolved Pho- 
toemission Spectroscopy (ARPES) nor the evolution of 
the transport phenomena across the HTC phase diagram. 
The MFL hypothesis assumes that optimum T c lies in 
proximity to a quantum critical point and consequently, 
quasi-particles are ill defined everywhere on the Fermi 
surface with a self-energy that is governed simply by 
the temperature scale. In contrast to the two-lifetime 
picture, this model is consistent with certain aspects of 



ARPES measurements, though not all. Most notably, 
the recently reported evidence of a kink in the energy 
dispersionE3 seems at odds with a scenario based on quan- 
tum criticality. Yet again, predictions of the MFL theory 
have so far been confined to optimum doping. 

In this paper, I propose a new model for the normal 
state scattering rate that offers an alternative route to 
understanding the HTC phenomenology. The salient fea- 
tures of the model are presented in section II. The model 
has as its basis, two novel but physically reasonable as- 
sumptions. Firstly, the dominant scattering rate is as- 
sumed to vary as T 2 everywhere on the Fermi surface, 
but in contrast to previous approaches, both the elastic 
and the inelastic parts of the self-energy contain strong 
(doping dependent) basal plane anisotropy. The second 
key assumption is that this anisotropic scattering rate 
eventually saturates (on different regions of the Fermi 
surface at different temperatures) at a value equivalent 
to the Mott-Ioffc-Regel criterion for coherent band trans- 
port. In section III A, the model is compared with ex- 
perimental data on Tl2Ba2CuOg+5 (T12201) at different 
doping concentrations. Excellent qualitative and quanti- 
tative agreement is found between the model and exper- 
imental data for p a b(T), cotO#(T) and significantly, the 
magnitude and T-dependence of Ap a b/ p a b(T), with rel- 
atively few fitting parameters. The model goes beyond 
previous proposals in that it can also account for the evo- 
lution of the observed normal state behavior across the 
overdoped (OD) region of the phase diagram from op- 
timum doping. In section III B, the model is applied 
to recent ARPES measurements of the quasi-particle 
scattering rate and normal state transport data for 
I^S^CaC^Os+a (Bi2212) and comparisons are made 
with the single-layer Bi2Sr2~ :E La a ;Cu06+ ( 5 (Bi2201). In 
particular, the model is shown to provide a natural ex- 
planation for the non-universal power law -behavior of 
cot8 h(T) observed in the Bi-based cuprates .E3Ej Finally, 
in section IV, justification arguments for this particular 
choice of scattering rate are presented and the implica- 
tions of this phenomenological model for the physics of 
HTC across the phase diagram are discussed, with sug- 
gestions for future study. In the process, I challenge two 
of the previously established paradigms of the normal 
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state in HTC; namely the linear low-T behavior of p a b(T) 
and the relevance of the non-saturating p a b(T) observed 
at higher temperatures. 

II. MODEL 

There_is strong experimental evidence from both 
ARPESlia and c-axis MIO that the dominant scattering 
mechanism in HTC is highly anisotropic within the basal 
plane, being strongest at the so-called 'anti-nodal' points 
at (tt, 0) and weakest at the 'nodal' points (tt, tt). In an 
anisotropic 2D metal, Ong has shown that the Hall an- 
gle tan0#(T) will be dominated by those regioas of the 
Fermi surface where the scattering is weakest 113 If the 
anisotropy is sufficiently high, cotOir(T) should there- 
fore reflect the T-dependence of the scattering rate at 
the nodal points. Given that in the majority of cuprates, 
cot9 h(T) ~ A + BT 2 , I take this to represent the in- 
trinsic T-dependence of the dominant scattering rate in 
HTC. A key stipulation of the present model is that the 
scattering rate has this same T-dependence everywhere 
on the Fermi surface, but rhi*Sn four-fold symmetry, be- 
ing maximum along (n, 0).t3O Combined with an elas- 
tic scattering rate To((j)), which may also contain 4- fold 
anisotropy,E2l the 'ideal' scattering rate Pideai(0) is writ- 
ten as 

r idea i(</»,T) =r o (0) + r r (^T) 

where 

r (» = a(l + ccos 2 2(f) 
and the T-dependent scattering rate 

T T (>, T) = 0(1 + e cos 4 2<j))T 2 

Here <p represents the angle between the in-plane Fermi 
wave vector and (k x , 0) whilst c and e are the anisotropy 
factors for the impurity and T-dependent scattering rates 
respectively. The coefficients a and (3 can be roughly es- 
timated from Hall angle measurements. The particular 
angular dependencies for Po(</>) and Tt(4>) shown above 
are chosen to reflect the four- fold symmetry of rid e ai(</>) 
and to be consistent with ARPES measurements, as out- 
lined in section III B. Pideai(^) is largest (= a(l + c) + 
(3(1 + e)T 2 ) along (tt, 0) and smallest, but still finite (= 
a + (3T 2 ), along (tt, tt). 

The second key feature of the model is the application 
of the Mott-Ioffe-Regel (MIR) limit to HTC. In its ba- 
sic form, the MIR criterion states that the quasi-particle 
mean free path £ cannot become smaller than the lat- 
tice spacing a, since at this point, the quasi-particles lose 
their coherence and conventional Boltzmann transport 
analysis becomes irrelevant. When fc-dependence is in- 
troduced, different regions of the Fermi surface become 
incoherent at different temperatures. Likewise, the scat- 
tering rate at different points on the Fermi surface will 
saturate (at the MIR limit) at different temperatures. 



The first evidence for resistivity saturation in metals, 
consistent with the MIR limit, was reported in the A- 15 
superconductors.EEl The T-dependence of the resistivity 
was found to ftL. extremely well to a simple 'parallel- 
resistor' model,Efl suggesting that the ideal resistivity 
(i.e. in the absence of saturation) was somehow shunted 
by a large saturation resistivity p max corresponding to 
£ = a. Using, the MIR criterion as his guide, Gurvitch 
later arguedcj that there must be a minimum time be- 
low which no scattering event can take place, i.e. a mini- 
mum distance a over which each carrier gains additional 
drift velocity before being scattered, and from this de- 
rived an elegant understanding of the appropriateness of 
the parallel-resistor model for resistivity saturation. 

In the spirit of Gurvitch's picture, I introduce for HTC 
a maximum scattering rate (minimum lifetime) r max (</>), 
equivalent to £ — a, and define an effective scattering 
rate r c ff (0, T) at each point on the Fermi surface of the 
parallel-resistor form, i.e. 

1 _ 1 1 

r e fi(0, t) r idca i(0, t) r max (^) 

T g((f>, T) is then inserted into the Jones-Zener ex- 
pansion of the standard Boltzmann transport equation 
for a quasi-2D Fermi surface to extract all measurable 
transport coefficients. A derivation of the appropriate 
expressions is included as an Appendix. Since the Fermi 
velocity vp (= 1 - 3 x 10 5 ms _1 ) has already been deter- 
mined by ARPES, the value of r max can be fixed directly 
by the MIR criterion (£ = a) giving r max in the range 
250 - 750 THz (150 - 450 meV). As will be shown, this 
derivation of r c ff(0, T) is sufficient to explain the many 
unusual T-dependencies that appear in the normal state 
transport in HTC. 

The T-dependencies of T e g((f>, T) for the two momen- 
tum directions (tt, 0) and (tt, tt) are illustrated schemat- 
ically in Fig. |l]. For concreteness, the parameters used 
in this simulation are those appropriate for OP T12201, 
namely a = 18 meV, (3 = 1.15 x 10~ 3 meV/K 2 , c = 6, e 
= 30 and r ma >'°) = 435 meV. Along (tt, 0), r id cai(T) 
reaches r max at T c ~ 90K, at which point, the anti- 
nodal quasi-particles will become incoherent. Due to the 
presence of the 'shunt' r max , however, r o ff(T) is always 
smaller than Fjdeai(T) and approaches saturation much 
more gradually. For the quasi-particles at (tt, tt) on the 
other hand, Fidcai(T) reaches r max at a much higher tem- 
perature (T ~ 600K). Hence, at the nodal points, well- 
defined coherent quasi-particles exist at all relevant tem- 
peratures. Note that r e ff(T) at (tt, tt) is quasi-linear over 
a wide temperature range. The 'effective' anisotropy fac- 
tor e* (= reff^' ) /r off - 1) is significantly lower than 
e for all finite T, reaching a maximum value ~ 7 at the 
temperature where Tideai^' '' first crosses r max (see inset 
to Fig. 111). The overall effect is a transport scattering 
rate T e g(cp, T) whose anisotropy initially grows with in- 
creasing T but then gradually becomes more isotropic at 
higher temperatures as the scattering rate at different re- 
gions of the Fermi surface tends towards saturation. In 
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FIG. 1. T-dependence of the ideal and effective scattering 
rates for (tt, 0) and (n, n) within the model. The parame- 
ters used are given in the text. Inset: T-dependence of the 
effective anisotropy e* = rW (,r ' 0) /r o ff ( ' r '' r) - 1 for the same 
parameters 

this way, the T-dependence of e*(T) mimics that of the 
Hall coefficient R H (T) in HTC. 



III. FITTING 

A. Application to T12201 

T12201 is an ideal system for modelling the normal 
state transport properties in HTC for several important 
reasons; (i) T12201 is a single band cuprate, and there- 
fore does not contain any secondary conducting subsys- 
tem, such as the CuO chains that complicate Hall and/or 
MR data in YBa 2 Cu 3 07_5 (YBCO); (ii) having a sin- 
gle Cu02 plane, bilayer coupling is irrelevant; (iii) the 
doping range of T12201 covers only the OD region of 
the phase diagram thereby avoiding any further compli- 
cations due to the normal state pseudogap; (iv) a T- 
dependent anisotropy in the in-plane scattering rate has 
already been observed in £p T12201 (T c = 30K) from 
c-axis MR measurements ;E3 and finally (v) p a b(T) and 
cot©#(T) exhibit the very clean linear and quadratic T- 
dependencies characteristic of OP cuprates as shown in 
Fig. U (±he_experimental data are represented by closed 
circles) t3t3 One surprising feature of the experimen- 
tal data is the small negative zero-temperature intercept 



in Pab(T) (extrapolation indicated by a dotted line) co- 
existing with a large positive intercept in cotO h(T). Line 
fits to the experimental data between 130K and 300K ac- 
tually give p ab (T) = -10 + 1.56Tluflcm) and cot9 H (T) 
= 27 + 0.0014T 2 (in 10 Tcsla).B This subtle but im- 
portant feature has often been overlooked but as will be 
shown, is a natural consequence of the form of the scat- 
tering rate presented here. 

To complete the Boltzmann transport analysis, I in- 
troduce a Fermi wave vector kp and Fermi— velocity v F 
consistent with band structure calculations ,E2l i.e. fcp(</>) 
= 6.5(1 + 0.15 sin 2 2</>) A" 1 and v F (4>) = 2.5(1 + 0.2 
sin 2 20) x 10 5 ms _1 . This expression for kp((f>) is consis- 
tent with a hole doping concentration p — 0.16. The small 
anisotropy factors do not play a significant role in the fol- 
lowing analysis and for simplicity, both kp(cj)) and vp((f>) 
are assumed to be independent of doping concentration. 
The corresponding saturation scattering rate r max (</>) = 
650(1 + 0.2 sin 2 20) THz or 435(1 + 0.2 sin 2 20) meV 
(since a = 3.86A). Note that the value of r max is en- 
tirely determined by our choice of vf and contains the 
same anisotropy factor as vp to ensure that the condi- 
tion I = a is satisfied everywhere. Again, these values are 
kept constant for all doping concentrations investigated 
here. Thus in the ensuing analysis, there are only four 
adjustable parameters, namely the two anisotropy pa- 
rameters c and e and the coefficients of Fjdeai^), a and 
(3. The choice of a and (3 is somewhat restricted how- 
ever since according to the Ong representation, a and 
P should be comparable with the coefficients A and B 
determined from cotOj^(T). These four parameters are 
then adjusted to obtain the best global fits, labelled here- 
after p${T), cote^(T) and Ap f a f /p f a f(T), to the three 
independent transport properties p a b(T), cot@#(T) and 
Ap a b/ p a b(T). Finally, scaling corrections of the order 
30% were allowed in the fits to account for any system- 
atic errors in the experimentally determined quantities. 

The solid lines in Fig. | are p f a f(T) and cot6^'(T) 
for OP T12201. The fitting parameters used are a = 27 
THz, j3 = 0.00175 THz/K 2 , c = 6 and e = 30. Above 
130K, both the linear resistivity and the quadratic Hall 
angle are reproduced with slopes and (extrapolated) T 
= intercepts in excellent agreement with the experi- 
mental data. Below 130K, p„j'(T) deviates from linear- 
ity, approaching T 2 as T — > 0, whilst p a b{T) actually 
shows a downward deviation from linearity below 130K, 
presumably due to the onset of superconducting fluctu- 
ations which always dominate the intrinsic behavior of 
the normal state p a b(T) near T c . Note that the T depen- 
dence of r idcal (T) (= 27 + 0.00175T 2 ) is similar to that 
of cotO H (T) (p. 27 + 0.0014T 2 ), in good agreement with 
Ong's pictureO 

The insets in Fig. |2| show p^f(T) and cotO^*(T) over 
an extended temperature range up to 500K. The influ- 
ence of r max on the high-T behavior is manifest as a slight 
deviation from linearity above 360K, indicating the on- 
set of saturation in r o ff(T). The fact that this tendency 
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FIG. 2. Topf-p a( ,(r) (closed circles) and p£"(T) (solid line) 
for OP T12201.B The dotted line is an extrapolation of p a b{T) 
to T = OK. Bottom: ce*@ H (T) (closed circles) and cot9^"(T) 
for the same crystal.Eil Insets: The same fits up to 500K, 
showing the onset of saturation at higher T. 



towards saturation is not actually observed in p ab -(T) in- 
dicates that the model does not capture all the essential 
physics. However, as I argue in section IV, there are sev- 
eral possible factors that might lead to a non-saturating 
resistivity in HTC which go beyond conventional Boltz- 
mann transport analysis and may mask the intrinsic be- 
havior of the scattering rate at higher temperatures. 

The modified Kohler's rule Ap a b/p a b = m 2 tan 2 0# 
is also— obeyed in OP T12201, with m ~ 1.9 down to 
130K,E3 as shown by the dashed line in Fig. |^. (Be- 
low this temperature, superconducting fluctuations again 
begin to dominate.) This value lies in between the ex- 
perimentally determined .values for OP La2- :c Sr :z: Cu04 
(LSCO)0 and OP YBCOB indicated by dotted lines in 
Fig. | The open circles in the figure are m' lt (T) for 
OP T12201 calculated using the same model parameters 
that were used in Fig. |^. As can be seen, mf lt (T) has 
the correct magnitude and T-dependence, decreasing by 
only 10% between OK and 300K. Note that the present 
model does not include vertex corrections nor paracon- 
ductivity contributions that might influence the form of 
Apab/Pab{T) at low T. 

The behavior of m^ lt {T) is determined almost exclu- 
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FIG. 3. Dashed line: Modified Kohler'srfdot m(T) for OP 
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ybco.eI 



sively by the anisotropy parameters c and e. (If c = 0, 
for example, Ap^f /plf(T) would fall steeply as T — > 0.) 
The large value of e is surprising, but is necessary, not 
only to give the correct magnitude for m at finite tem- 
peratures, but also to maintain a robust T 2 dependence 
of cotO^*(T). Hence, within this model, the observation 

of a quadratic Hall angle and the form of Ap^ / p^f (T) 
are all signatures of an extremely anisotropic basal plane 
scattering rate that likewise varies as T 2 . The T-linear 
resistivity, on the other hand, is simply an artefact of the 
broad crossover between the (approximately) T 2 resistiv- 
ity at low T and the saturation of r c ff at high T. 

The ability of the present model to reproduce the cor- 
rect magnitude and T-dependence of the orbital MR is 
to be contrasted with the 'Cold Spots' modelQ where the 
parameters required to separate the transport and Hall 
lifetimes yield an orbital MR that is orders of magni- 
tude too large and has a much stronger T-dependence 
than shown in Fig. ^. This difficulty is overcome in the 
present model by the introduction of the shunt r max that 
makes the effective anisotropy around the Fermi surface 
much smaller (see inset to Fig. |j]), thereby reducing both 
the magnitude and the T-dependence of Ap^f /p^f(T). 

The evolution of the normal state transport properties 
in T12201 with doping are summarized in Fig. |]. The top 
panel in Fig. § shows p a b{T) (closed circles) and p^iT) 
(solid lines) for two OEL.T12201 single crystals with T c 
values of 55K and 30Kj£3 The parameters used in these 
fits are a = 9 THz, (3 = 0.0016T 2 THz/K 2 , c = 2 and e = 
11.5 for the 55K crystal, and a = 7 THz, (3 = 0.0015T 2 
THz/K 2 , c = 2 and e = 8 for the 30K crystal. Note that 
the anisotropy factors are markedly reduced compared to 
OP T12201, reflecting the increased curvature observed 
in p a b(T). The bottom panel in Fig. ^ shows the corre- 
sponding cotO^f(T) and cotO^*(T) for the two samples. 
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FIG. 4. Top: p ab (r)Ei (closed circles) and (solid 
lines) for,«D T12201 (T c = 55K and T c = 30K). Bottom: 
cotO h(T)E3 (closed circles) and cotO^ !t (T) (solid lines) for 
the same samples with the same fitting parameters. Inset: 
m^ lt (T) (open circles) for the two sets of parameters. The 
dashedJine represents the experimental data for the 30K 
crystal £3 



In contrast to OP T12201, cot9§*(T) deviates slightly 
from a simple T 2 dependence. In fact, the T-dependence 
is closer to T 1 - 95 and T 19 for the 55K and 30K crys- 
tals respectively. This deviation from a T 2 dependence 
reflects the fact that in less anisotropic cuprates, other 
regions of the Fermi surface are beginning to contribute 
to the Hall angle. The inset in the bottom panel of Fig. |] 
gives the modified Kohler's plot derived for the two dop- 
ing levels. As we expect, m^ lt is lower in these examples 
compared with OP T12201, with values approaching 1.55 
and 1.4 for the 55K and 30K crystal respectively. Experi- 
mentally, to is found, to be 1.4 for 30K T12201 (indicated 
by a dashed line) ^3 again in excellent agreement with 
the model. 



B. Application to Bi2212 

With its improved energy and momentum resolution, 
ARPES has become a powerful probe of the single- 
particle self-energy S of HTC. The vast majority of 
work to date has been carried out on the bilayer cuprate 
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FIG. 5. Quasi-particle scattering pate for OP Bi2212 ex- 
tracted from ARPES measurements Ej The various numbers 
and symbols refer to different locations on the Fermi surface 
as described in the text. The solid lines are fits to the indi- 
vidual scattering rates using the expression for T e g(<j>, T). 



Bi2212, due to its ease of cleaving. One of the most strik- 
ing features of the emerging spectral function in Bi2212 
is the broad featureless normal state spectrum at (tt, 
0) whose lineshape sharpens dramatically into the well- 
known peak-dip- hump structure upon cooling below T c . 
Although interpretation of ARPES spectra is still a con- 
troversial topic,E-3 it has been arguecocJ that this dra- 
matic rearrangement of the ARPES lineshape indicates 
the presence of ill-defined quasi-particle states in the nor- 
mal state at (tt, 0), rendered incoherent by some form of 
catastrophic low-energy scattering mechanism. ARPES 
data also reveal that the broadening of the quasi-particle 
self-energy is much less dramatic for quasi-particles at (tt, 
tt) and that this basal plane anisotropic is much larger for 
OP cuprates than for OD cuprates.t3 All these features 
indicate a highly anisotropic normal state scattering rate 
whose anisotropy grows weaker with overdoping, consis- 
tent with the picture described above for T12201. 

Fig. H shows the variation of the quasi-particle scatter- 
ing rate with angle and temperature for an OP Bi2212 
single crystal, extracted from the ARPES data of Valla et 
al.E3 The various symbols shown in Fig. ^ are the same 
as those used in the original Valla plot of the momen- 
tum widths Afc, obtained in turn by fitting different mo- 
mentum distribution curves (MDC) with Lorentzian line 
shapes. The quasi-particle scattering rates 1/r are ob- 
tained using the usual expression 1/r = |2ImS| = HAkvk 
where ImE is the imaginary part of the single-particle self 
energy and Vk is the unrenormalized band velocity ~ 4.0 
x 10 5 ms _1 . The numbers refer to locations on the Fermi 
surface at angles, (1) = 45°, (2) = 37.5°, (3) = 25.5°, (4) 
= 17.5° and (5) = 7° away from (tt, 0). 
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FIG. 6. Angular dependence of To (top panel) and Yt (bot- 
tom panel) for OP Bi2212 determined from the fits to the 
MDC curves shown in Fig. |j| The solid lines are best fits 
whose expressions are shown in the appropriate figures. 



The normal state Fermi velocity, determined at each 
location from the slope of the dispersion relation near 
the Fermi level, can be approximated by the expression 
V F (<p) = 2.5(1 + 0.2 sin 2 20) x 10 5 ms-'O Note that 
this is the same vf(4>) used in the T12201 analysis. Simi- 
larly, r max (</>) = 650(1 + 0.2 sin 2 2</>) THz or 435(1 + 0.2 
sin 2 20) meV. 

The solid lines in Fig. | are the fits to the ARPES 
data using the model scattering rate r e ff(</>, T). The fits 
are all reasonable and give angular dependencies for the 
resulting Po(0) and Tt(<P) as shown in Fig. |(| The solid 
lines in Fig. g are weighted fits to the expressions To(<f>) 
= a(l + ccos 2 (20) and T T (<f>) = (3(1 + ecos 4 (2<?(>) with 
a = 140 THz, p = 0.0065 THz/K 2 , c = 4 ± 1 and e = 
8 ± 1. Note that Tt(4>) 1S more strongly peaked at (n, 
0) than Tq ((/>). The reason these anisotropy factors seem 
large is the presence of r max . As T e g(cf>, T) approaches 
Tmax, rideai(</', T) must change significantly in order to 
create even a small change in r o ff(0, T). Hence within 
this model, the 'ideal' scattering rate is actually much 
more anisotropic than the raw data suggest. 

The form of T c jf((p,T) obtained above is inserted into 
the Boltzmann equation to obtain a crude estimate of 
f lt (T) and cotO^*(T) as shown in the top panel of Fig. 



The resistivity has a very large zero-temperature offset 
though is approximately linear from T c up to T = 250 K 
at which point, pJuiT) tends towards saturation as the 
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FIG. 7. Top panel: p a b(T) and cotOff(T) (inset) for OP 
Bi2212 using the fitting parameters derived from the ARPES 
data. Solid circles: the results from the model, dashed lines: 
guides to the eye. Bottom panel: Best fits to p a b(T) (main 
figure) cote H (T) (top inset) and l/(Ap ab /p ab )°- B (T) (bottom 
inset) for OP Bi2212. Note the different T-dependencies of 
the various plots. 



effect of r max becomes prominent. cotO^*(T), on the 
other hand, is found to vary roughly as A + BT 15 . Ex- 
perimentally, p a b{T) in OP Bi2212 shows the ubiquitous 
T-linear behavior, but has a positive zero-temperature 



intercept, indicative of increased inhomogeneity in this 
compound,E3-flthilst cote^(T) = A + BT n where n = 
1.75 ± 0.05, EM in contrast to the usual result for OP 
cuprates-jii. = 2. Similar behavior is also observed in 
Bi220lMiy 

Whilst the ARPES-derived fits are in reasonable agree- 
ment with the experimental data, they are not ideal. In 
any case, one does not expect the transport scattering 
rate to be equivalent to the single-particle scattering rate 
obtained by ARPES, since the former should be domi- 
nated by large-angle scattering events whilst the latter 
will not. The best actual fits to the transport data for 
OP Bi2212E3 are shown in the bottom panel of Fig. ^ for 
the fitting parameters a = 50 THz, /3 = 0.00325 THz/K 2 , 
c = 4 and e = 8. These values give the expected linear 
resistivity between T c and 300K with the correct slope 



and intercept and cot6£ (T) = A + BT 1 ' 5 over the 
entire temperature range, again with the correct slope 
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and intercept. The deviation of cotO#(T) from its sim- 
ple T 2 dependence arises from the increased disorder in 
Bi2212 that must somehow smear out the anisotropy of 
the scattering rate around the Fermi surface. Note that 
the anisotropy factors are exactly the same as those ob- 
tained from the ARPES data, while (3 in the transport fit 
is half the ARPES value. The much larger residual scat- 
tering rate inferred from the ARPES spectra is a common 
feature of photoemission dataliJ so is not unexpected. In- 
cidentally, using the same parameters, Ap^ / p^£ (T) is 
found to vary as (C + DT 19 )- 2 , also shown in Fig. @. 
This is only a prediction from the model at present since 
Apaf,/ p a b(T) has yet to be established experimentally in 
OP Bi2212. 

It should be stressed .-finally that the ARPES mea- 
surements of Valla et alSn are often cited as strong evi- 
dence for a linear quasi-particle scattering rate in HTC 
that is itself responsible for the T-linear resistivity ob- 
served at optimal doping. However, as noted by Valla 
and co-workers, near the anti-nodal point (label num- 
ber (5) in Fig. ||), the scattering rate becomes almost 
T- independent. This tendency towards saturation is an 
essential part of the present model, and combined with 
a single scattering rate that is quadratic in temperature 
everywhere on the Fermi surface, it is able to account 
for the quasi- linear scattering rate at the nodes (Fig. ||), 
the T-linear resistivity, the tendency towards saturation 
at (tt, 0) and an inverse Hall angle that varies as A + 
BT 1 - 75 , without requiring any additional assumptions. 
Note that in the Cold-Spots model, the scattering rate 
at the nodes is always quadratic, whilst in the MFL phe- 
nomenology, there is no saturation at (tt, 0). 

IV. DISCUSSION 

It is clear that the physical mechanism behind the 
normal state scattering rate, and ipso facto high- 
temperature superconductivity, can only be determined 
once a coherent description of all the normal state trans- 
port properties in HTC is found. In this paper, I have 
outlined a new phenomenology that shows promising 
signs of achieving this goal. However, whilst many of 
the anomalous features can be explained within this sim- 
ple model, the current proposal appears, at first sight 
at least, to be incompatible with two key aspects of the 
experimental situation in HTC; namely the linear resis- 
tivity at low T and the absence of resistivity saturation at 
high T. Therefore, before going on to discuss the possible 
implications of the model, it is perhaps worth outlining 
to the reader how these two major discrepancies might 
possibly be reconciled. 

It has often been claimed that p a b(T) in OP cuprates 
is linear at all finite temperatures. This long- held view 
is based largely on the observation back in 1990 of a lin- 
ear p a b{T) in single-crystal Bi2201 extending from 700K 
down to T c = 10K.E3 If this were indeed the intrinsic be- 



havior, it would certainly require an explanation outside 
the realms of conventional FL physics. Important new 
data on the latest generation of Bi2201 crystals (with 
significantly improved T c values) tell a rather different 
story however. By suppressing superconductivity in OP 
Bi2201 (T~j= 32K) in a 60 Tesla pulsed magnetic field, 
Ono et alS3 recently found that the linear p a b(T) seen 
at high T (and giving an apparent zero intercept at T = 
OK) crosses over to a higher power T-dependence below 
around 60K, before saturating at some finite positive p*u 
Similar behavior has also been reported in OP LSCO.E3 
The extension of the linear p a b(T) to low temperatures 
in Bi2201 and LSCO (in zero-field) is not therefore in- 
trinsic and is presumably caused by the onset of fluctu- 
ation conductivity, as is evident in the p a b(T) data of 
T12201 and Bi2212. In addition, the (extrapolated) neg- 
ative intercept in p a LOJX found in the best quality T12201 
and YBCO crystalaSJa is unphysical and dictates that 
Pab(T) must cross over to a higher power of T at lower 
temperatures. Such behavior is incompatible with the 
MFL hypothesis, for example, but is a natural conse- 
quence of the form of T c g((j), T) presented here. Perhaps 
in the light of these new results, the long-standing as- 
sertion that p a b{T) obeys a simple linear T-dependence 
at low T should now be re-examined. It would certainly 
be interesting to follow the T-dependence of OP cuprates 
with higher T c to lower temperatures, as and when larger 
magnetic fields become available. 

I now turn to address the second point, namely the 
absence of resistivity saturation at higher temperatures. 
The MIR limit and the parallel-resistor model have been 
applied successfully to a wide range of materials over the 
years. Their relevance to strongly correlated electron sys- 
tems such as HTC, however, has been challenged in re- 
cent timcs.E-3 In these so-called 'bad metals', p(T) does 
not saturate. Instead p(T) grows approximately linearly 
with T at high temperatures (typically T > 300K) reach- 
ing values of the order 1-10 mf2cm that correspond 
to i <C a. It is claimed that this absence of resistiv- 
ity saturation is a signature of a non-FL ground state 
at zero tcmperature.E-a At first sight, these arguments 
appear rather compelling. However, the origin of the 
non-saturating resistivity in these systems is still poorly 
understood. More importantly, similar p(T) behavior has 
also been reported in S^RuO^lS a close structural ana- 
log of the cuprates which exhibits quantum oscillationsE3 
and so forms a well-defined FL ground state at low T. 
The association of non-saturating p(T) with a non-FL 
ground state is therefore misleading and that some other 
physical origin must be sought for the high-T behavior. 

Possible sources for the excess resistivity include ther- 
mal expansion effects (known to play a significant role 
in the organic conductors, for example), a thermally- 
induced reduction in the density of states or carrier con- 
centration, and the onset of a 'thermal' diffusive (non- 
Boltzmann) transport regime at high T that sets in once 
the concept of a scattering rate is no longer applica- 
ble and temperature becomes the only relevant energy 
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scale in the problem. Determining the origin of this 
non-saturating resistivity, though desirable, is not essen- 
tial for justifying the model however, since the current 
model is only strictly applicable in the low-T (T < 300K) 
regime, where the majority of the carriers are still well- 
defined quasi-particles and Boltzmann transport analysis 
remains valid. 

Any observation of a threshold in the scattering rate 
would obviously lend significant weight to the arguments 
presented here however. In addition to the work of Valla 
et aL,E3 Norman et al. have recently extracted ImS 
along (tt, 0) for slightly OD Bi2212E3 and found that 
it too saturates above T c at a value 2Im£ max ~ 180 
meV, well within the range of validity given by I = a. 
Optical measurements can also be used in principle to 
separate T from the conductivity, though here the situa- 
tion is not as clear. In particular, T is found to saturate 
with increasing temperature, but NOT with increasing 
frequency.EZI Finally, recent analysis of the plasmon spec- 
tra in K3C60J23 another system with non-saturating re- 
sistivity, has revealed that the quasi-particle scattering 
rate itself saturates around T — 500K at a value r max ~ 
450 meV. Whilst more experimental work is clearly re- 
quired to clarify this effect, these observations do point 
to the possibility that the MIR limit may still be a pre- 
serve of all metallic systems, in the strictest sense that 
r never exceeds the maximum value r max corresponding 
to £ — a, and that the excess resistivity seen in these 
strongly correlated systems should not be regarded sim- 
ply as a signature of an ever-increasing scattering rate. 
Moreover, since the phenomenology described in this pa- 
per has been shown to provide a coherent description of 
many of the anomalous normal state transport properties 
in HTC, further efforts to separate T(T) from p(T) are 
strongly recommended. 

Turning now to the implications of the model, the most 
striking features of the analysis on T12201 are clearly 
the large values for c and e and their strong variation 
with doping, as summarized in Fig. |8|. (For each crystal, 
the doping level p is obtained from the well-knpwn HTC 
parabola law T c (p)/T c max = 1 - 82(p - 0.16) 2 ,El shown as 
a solid line in Fig. ||.) 

The large value of c, particularly at optimum dop- 
ing, is somewhat surprising, though it is not inconsistent 
with the picture emerging from ARPES-jpeasurements 
on Bi2212B3 Recently it has been arguedEl that scatter- 
ing on impurities located off the Cu02 planes may be 
the source of this four- fold anisotropy in To(<fi). In the 
case of tetragonal T12201, such interplane randomness 
might be induced by Tl/Cu site substitutions as well as 
by interstitial oxygen. It is not clear at this stage why 
such large anisotropy would also appear in the transport 
scattering rate, which typically depends only on large- 
momentum— transfers. However, as noted by Varma and 
Abrahams,tj such anisotropic elastic scattering may well 
play a key role in the magneto-transport behavior. In- 
deed, the large value of c is only required in the fitting 
here to maintain the constancy of m^ lt (T) down to low 
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p Doped Holes/Cu 

FIG. 8. Variation of the anisotropy factors c (open dia- 
monds) and e (closed diamonds) with doping concentration 
p. The dashed lines are guides to the-eye. The solid line rep- 
resents the T c (p) parabola for HTC.E3 The error bars reflect 
the range of fitting parameters that provide a reasonable fit 
to all three transport coefficients. 

temperatures. 

The T-dependent anisotropy factor e is found to drop 
steeply from optimal doping, scaling approximately with 
T c and approaching zero close to the doping level where 
superconductivity vanishes. This is to be expected per- 
haps, since the separation of transport and Hall lifetimes, 
is observed to vanish in the most highly OD samples.c£l 
The correlation between e(p) and T c {p) on the OD side is 
intriguing and suggests that the mechanism giving rise to 
this strongly anisotropic scattering rate might also be as- 
sociated with the anisotropy (d-wave pairing symmetry) 
of the superconducting order parameter. 

As illustrated in Fig. |, T idc J- nfi \T) for OP T12201 is 
found to cross the MIR limit close to T = T c . (Fe^* ) (T) 
will obviously be lower than the MIR limit due to the 
presence of the shunt r max .) The same is also true for 
OP Bi2212. The picture that emerges therefore is one 
in which the quasi-particles at (tt, 0) are approaching an 
incoherent state at T = T c , consistent with what is ob- 
served in ARPES, and this may well have implications for 
the onset temperature of superconductivity on the under- 
doped (UD) side of the phase diagram. More explicitly, 
if e continues to grow with decreasing p, the temperature 
at which the anti-nodal quasi-particles become incoher- 
ent will decrease accordingly. What ultimately causes 
the suppression of T c on the UD side of the phase di- 
agram is still an open issue, and I cannot address that 
important question here, since clearly, the physics of the 
UD cuprates is also heavily influenced by the opening of 
the pseudogap which cannot yet be accounted for within 
this simple phenomenological model. Nevertheless, the 
onset of incoherent quasi-particles at T = T c max is an in- 
triguing observation from the above analysis, hints at a 
new interpretation of the T c (j>) parabola across the HTC 
phase diagram and should be investigated further. 
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FIG. 9. Simulation of the variation of n(p) for Bi2212 
where n is the coefficient of inverse Hall angle cotO h(T) — 
A + BT n . The anisotropy factors (c, e) used in the fitting 
routines are indicated in parentheses. The dashed line is a 
guide to the eyeO Inset: Effective anisotropy e*(T) for two 
selected points. 

Some possible insight into the nature of the scatter- 
ing on the UD side can be gained from further analysis 
of the Hall angle in Bi2212. In recent detailed trans- 
port studies of the Bi-based cuprates Bi2212 and Bi2201, 
cotO/f(T) was found to vary as T n with n decreasing 
from n ~ 2 to n ~ 1.6 - 1.7 as one moves from the 
UD to the OD regime. This non-universal power law 
behavior in cotOff(T') reveals a new layer of complex- 
ity in the phenomenology of normal state transport in 
HTC that has yet to be addressed. It was noted earlier 
for T12201, that n decreases on the OD side as c and e 
are reduced. Given that at optimal doping, n ~ 1.8 lies 
roughly in the middle of the variable range for Bi2212, 
it is tempting to attribute this variation in n(jp) purely 
to changes in c and e across the phase diagram. A sim- 
ple example of how this might come about is illustrated 
in Fig. H where c and e have been varied systematically 
to simulate the experimentally observed variation of n 
with doping]!] whilst keeping a and (3 fixed at the val- 
ues extracted from the transport fits for OP Bi2212 in 
Fig. ^. Intriguingly, in the isotropic case c = e = 0, n 
~ 1.65, in agreement with the experimental data. As c 
and e are increased, n rises as expected, but only ap- 
proaches n = 2 for extremely large values of e. Whilst 
this example is not completely rigorous, it does suggest 
that with more careful analysis, the systematic variation 
of n with doping in Bi2212 and Bi2201, and its range of 
values, might be understood within this simple model. 
It also implies that the anisotropy of the T 2 scattering 



rate must increase (almost exponentially) with decreas- 
ing doping on the UD side, reaching values greater than 
100 for p ~ 0.05! Indeed, within the present model, it is 
the only way to account for the observation of cotOn (T) 
= A + BT 2 in UD Bi2212. However, as illustrated in 
the inset of Fig. 0, the effective anisotropy e*(T) above 
100K is almost indistinguishable for the OP and the most 
UD cases. Hence, whilst e may increase dramatically, the 
more experimentally accessible quantity e* does not. 

Comparison of the fitting parameters for OP Bi2212 
and OP T12201 reveal another important aspect of the 
modelling. The degree of anisotropy in OP Bi2212 (c 
= 4, e = 8) is markedly reduced compared with that 
in OP T12201 (c = 6, e = 30) whilst the coefficients of 
ridcai(T) (a = 50 (27) and j3 = 0.00325 (0.00175) for OP 
Bi2212 (OP T12201)) are increased. This implies that 
the strength of the scattering rate, averaged around the 
Fermi surface, is very similar in the two cases, but the 
higher degree of disorder and inhomogeneity in Bi2212, 
evidenced by the larger value of a (the level of impurity 
scattering at the nodal points) acts to smear out this 
strong anisotropy. 

Whilst I do not wish to speculate at this stage on 
the origin of T e s((j), T), there are still some impor- 
tant remarks that can be made from the above analy- 
sis. Firstly, the strong anisotropy and T 2 dependence of 
T e g((f), T) strongly suggests that the dominant scatter- 
ing rate in HTC is electronic in origin and that electron- 
phonon scattering plays only a minor role in the nor- 
mal state transport properties of HTC. The large in- 
crease in the quasi-particle lifetime below T c observed 
in conductivity measurementsBc3 is consistent with the 
disappearance of this catastrophic electronic scattering 
mechanism below T c . Secondly, the large value for e for 
OP T12201 (Fig. ||) implies the existence of a singular 
scattering mechanism, peaked at (w, 0), whose inten- 
sity grows rapidly with decreasing carrier concentration. 
Moreover, the striking growth of the intensity of the scat- 
tering at (tv, 0) as we move towards the antiferromagnetic 
(AFM) insulating state (Fig. ||) points to a scattering 
mechanism that is magnetic in origin, possibly due to 
some form of AFM spin fluctuations. Another poten- 
tial candidate though, particularly given the form of the 
anisotropy and its singular nature, is the collective reso- 
nance mode seen in neutron scattering measurements on 
several cuprates, including T12201. Intriguingly, whilst 
the resonance mode is resolution limited, i.e. strongly 
peaked at (tt, 0) in T12201,Elit is significantly broadened 
in fc-space in Bi2212jfj consistent with the conclusions 
derived above. It is noted with caution however, that 
this mode has so far only been unambiguously resolved 
below T c , so its relevance to scattering in the normal state 
has yet to be established. 



9 



V. CONCLUSIONS 



VI. ACKNOWLEDGEMENTS 



In summary, I have demonstrated how the seemingly 
anomalous dc transport properties of HTC may be ac- 
counted for consistently within a Boltzmann transport 
analysis with a single T 2 scattering rate that varies 
markedly around the basal plane. In particular, the 
separation of transport and Hall lifetimes, the modified 
Kohlcr's rule and the evolution of the transport proper- 
ties across the OD region of the phase diagram can all be 
reproduced with a single set of parameters. In addition, 
the different T-dependence of the Hall angle observed in 
the Bi-based cuprates is shown to be simply an artefact 
of the increased inhomogcneity of these materials; the 
T 2 scattering rate continues to be the dominant source 
of scattering in Bi2212 but its singular nature is smeared 
out by the disorder. The form of T e s ((f), T) is also consis- 
tent with several features of the single-particle spectral 
function identified by ARPES, including the strong basal 
plane anisotropy (both in the zero-temperature and the 
T-dependent scattering terms) and the tendency towards 
saturation along (jr, 0). Such a wide-ranging consistency 
with experiment is indeed encouraging and strongly sug- 
gests that a conventional Boltzmann transport approach 
might be more appropriate to describe the normal state 
of HTC (at least in the intermediate and low T ranges) 
than was previously imagined. 

The fact that the parent HTC compound at half- 
filling is an AFM Mott insulator rather than a metal has 
pointed clearly to the role of strong electron correlations 
in the square planar Cu02 lattice as a key element in the 
cuprate problem. The introduction of holes into a highly 
correlated 2D magnetic background is believed to give 
rise to a non-FL ground state at low doping, character- 
ized by the highly anomalous physics of the UD cuprates, 
and ultimately, unconventional high temperature super- 
conductivity. It has often been assumed moreover, that 
with continued doping, HTC eventually evolve into a con- 
ventional FL as the electron correlations become weaker 
and the system becomes more three-dimensional. In this 
paper, I have identified a form for the normal state scat- 
tering rate that although somewhat unconventional, is 
still rooted in FL physics, yet can explain the evolution 
of the normal state transport properties of OD HTC right 
up to optimum doping. Of course, I do not expect this 
simple model to remain valid on the UD side of the phase 
diagram, but if it shown to be applicable to the OD and 
OP cuprates, then the way in which this relatively sim- 
ple picture breaks down, as the pseudogap opens and 
electron correlations become even stronger, might yet re- 
veal the crucial missing link between the FL fixed point 
at large doping and the curious physics of HTC on the 
underdoped side. 
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VII. APPENDIX: TRANSPORT COEFFICIENTS 
FOR A QUASI-2D FERMI SURFACE 

The current response Ji to an applied electric E is 
given by 

J * = J ev i9kd 3 k = (TijEj 

where g k is the displacement from the equilibrium 
distribution (g k = f - /o) due to the applied fields, ct^ 
is the conductivity tensor and Vi is the Fermi velocity in 
the i direction. The total displacement can be summa- 
rized in the Bloch-Boltzmann transport equation, which 
under the relaxation time approximation, is given by 

eE.v k — + -[v k xB}— = - 9k T 

The various components of erj™^ can be obtained 
through the Jones-Zener expansion 

to arrive at: 

To calculate, for instance, the in-plane Hall conductiv- 
ity a^y , we simply apply B z and E y and calculate the 
response J x using the Jones-Zener term g^ . 

For a quasi-2D Fermi surface, there are various tech- 
niques one can adopt to simplify the calculation. Firstly, 
the element d 3 k is the product of an area parallel to the 
Fermi surface dS and a radial (in-plane) component dk r 
which is related to an incremental energy de by 

de = hvp.kp — hvp cos 7dfc r 

Here 7 is the angle between vf and dk r (within the 
plane). Since vf is always perpendicular to the Fermi 
surface, simple geometrical considerations give 

70/0= tan- 1 (log M0))] 
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which is zero for a circular Fermi surface. 
The element dS can be expanded into the cylindrical 
elements kpd<p and dk z to give 



tan 6* = Pxv 



RhB 



u xy 

u XX 



271 



ir/d 



k, 



o J-7r/d hvpcoaj 



for a quasi-2D Fermi liquid since (-dfo/ds) is a delta 
function (provided all scattering is quasi-elastic). It 
is this integral J dcj) that will contain all the in-plane 
anisotropy of the essential parameters in our calculation. 
With B || c, the cross-product [vk x B] d/dk becomes 
Bvpd/dk\\ where 

d cos 7 d 
dkn 



kF d<j) 



whilst for B || ab, making an angle ip with the a-axis 
[v k xB}^ k =v z [zxB]^- + [v xy xB}-^- 



i.e. 



[v k x B] 



dk 



Bvp sin(0 — 7 — ip) 



dk. 



Finally, it is noted that the tetragonal symmetry of 
T12201 allows the limits of integration to be reduced to 
and 7r/2 by including an extra factor of 4. 

Employing these simple techniques, one can now derive 
the individual terms of the in-plane and c-axis conduc- 
tivity tensors. In this paper however, only the in-plane 
magneto-transport properties are considered and their 
formulae are listed below in the isotropic case. For the 
anisotropic case, simply replace Vf, kF and T by vf(<I>), 
fcp(</>) and T(<f>), or more appropriately for the present 
model, by T cS ((j), T). 



r(0) 



4^ ( T )4 io IWy d<t> 



r( 2 ) 



i B 2 ,2ti\. r /2 v F 



Att^ d' J r 



(-r) 4 / ^cos(0- 7 ) 



d v F cos-/ d v F . 
^ { ^-^ ( - C ° S( ^ 7))}d0 

The measurable quantities p a b, tan0^ and Ap a b/p a b 
are obtained by inversion of the conductivity tensor: 



Pab 



1 

V XX 



Pab 



a (2) 
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